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Abstract 

We calculate potentials between a proton and a H° (hyperon with strangeness —2) through the equal-time Bethe- 
Salpeter wave function, employing quenched lattice QCD simulations with the plaquette gauge action and the Wilson 
quark action on (4.5 fm) 4 lattice at the lattice spacing a ~ 0.14 fm. The ud quark mass in our study corresponds to 
~ 0.37 and 0.51 GeV, while the s quark mass corresponds to the physical value of rriK- The central pH° potential 
has a strong (weak) repulsive core in the 1 S'o ( 3 Si) channel for r < 0.6 fm, while the potential has attractive well at 
medium and long distances (0.6 fm < r < 1.2 fm) in both channels. The sign of the pH° scattering length and its 
quark mass dependence indicate a net attraction in both channels at low energies. 
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1. Introduction 

Modern nuclcon-nuclcon (TV AT) potentials pQ give successful description of the TVTV scattering data at 
low energies and have been used for precision calculations in light nuclei [2]. Similarly, hyperon-nucleon 
{YN) and hyperon-hypcron (YY) potentials in the coordinate space are known to be quite useful and 
important for studying the properties of hypernuclei [3j as well as the equation of state for hyperonic matter 
in neutron stars [4]. However, there are still large experimental uncertainties in YN and YY interactions at 
present, because the short life-time of the hyperons makes the scattering experiments difficult. Accordingly, 
phcnomcnological YN and YY potentials are not well constrained from data even under some theoretical 
guides [5I6I7I8I9H0] , 

In such a situation, it would be desirable to analyse hyperon interactions on the basis of the first principle 
lattice QCD simulations. Studies along this line for nucleon interactions was initiated in Ref.|ll|. where the 
AW scattering lengths were extracted from the quenched simulations by using the Luscher's finite volume 
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method [32]. The same method was applied later to the (2+l)-flavor simulations with the mixed action [T3"] . 
The scattering length of AN system was first examined in |14| in quenched simulations on a small lattice box. 
Subsequently, AN and S7V scattering lengths were studied in (2+l)-flavor simulations with the mixed action 
1 1 5 j . Such low-energy scattering parameters calculated in lattice QCD would be valuable inputs to construct 
phcnomcnological YN and YY potentials to be used for studying hypernuclci and hyperonic matter. 

Recently, an alternative but closely related approach to [12] has been proposed to define the NN potential 
from lattice QCD [16117118] , Since the potential is not a direct physical observable, one cannot make quan- 
titative comparison of this lattice NN potential with the phenomenological NN potentials. However, they 
are both designed to reproduce the correct scattering phase shifts, and their spatial structures are found to 
have common features [16] : the attraction at long and intermediate distances and the strong repulsion at 
short distance in the S'-wave channel |19j . 

The purpose of the present paper is to report our first attempt to apply the above approach to the hyperon- 
nuclcon systems. For the NN potential, only two channels (isovector and isoscalar) exist in the flavor SU(2) 
space, i.e., 2 <g> 2 = 3 © 1. Including the strange quark extends this to 8 ® 8 = 27 © 10* 0108© 10 ©8 
in the flavor SU(3) space. The isovector (isoscalar) channel in the NN sector is assigned to be a subset in 
the 27-plet (10*-plet) representation. The potentials in newly arising channels are hardly determined from 
experiments so far. 

In this paper, we focus on the iVS potential in the isovector (/ = 1) channel as a first step. (A preliminary 
account of this system has been reported in [20].) There are two main reasons for picking up this channel: 
(i) Theoretically, it is the simplest generalization of the NN system; pEP is obtained from pn by replacing 
the (i-quarks in the neutron by the s-quarks. Also it is a channel which do not have strong decay into 
other YN systems^] (ii) Experimentally, not much information has been available on the iVS interaction 
except for a few studies; a recent report gives the upper limit of elastic and inelastic cross sections |21j . 
and earlier publications suggest weakly attractive S-nucleus potential [22j . The S- nucleus interaction will be 
soon studied as one of the day-one experiments at J-PARC [23] via (K~ , K + ) reaction with nuclear target. 

This paper is organized as follows. In Section [2] we describe the basic formulation to derive the pS° 
potential through the Bcthc-Salpeter amplitude measured in the lattice QCD simulations. Our lattice setup 
is explained in Section [3] In Section [5] we show numerical results of the potentials in the spin-singlet and 
spin-triplet channels and their quark mass dependence. We also show the estimate of the scattering lengths 
in these channels. Section [5] is devoted to summary and concluding remarks. 



2. Basic formulation 



Our methodology to obtain the baryon potentials is along the lines of [16117118] . (See also [12124] for the 
seminal attempts to introduce similar notion of the potential.) We consider the low-energy iV-S scattering 
and start with an effective Schrodinger equation for the equal-time Bcthe-Salpetcr (BS) wave function 4>{r) 
obtained from the Lippmann-Schwinger equation [18|25| : 

1 V 2 0(r)+ / U(r,r')(f>(r')d 3 r' = E(f>(r). (1) 



2/' 

Here \i = m^ms/ '(mjy + m=) and E = fc 2 /(2/i) are the reduced mass of the JVS system and the non- 
relativistic energy in the center-of-mass frame, respectively Eq. (fTJ is derived from QCD by adopting local 
three-quark operator as the nucleon interpolating operator. It can be shown that the non-local potential 
U(r, f) is energy independent. Also, it is unique given the interpolating operator as long as the BS-amplitudc 
for all energies is known [18] . If we take other nucleon interpolating operator, both the BS wave function and 
the non-local potential change their spatial structure without affecting the obscrvablcs such as the binding 
energy and the phase shift. At low energies, the nonlocality of the potential can be expanded as [55] 

U(r,r , ) = V NS (f,V)5(r-r*), (2) 



Note that iVE in the isoscalar (/ = 0) channel is above the AA threshold. 
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where the dimensionless expansion parameter is V /M with M being the typical scale of the strong interaction 
such as the excitation energy of a single baryon. The general expression of the potential Vnb. is known to 
be [27] 

Vjvh = V (r) + V a {r){3 N ■ a s ) + V T {r)(T N • f H ) + V aT (r)(a N ■ ue.){t n ■ f H ) 

+V T (r)S 12 + V TT {r)S 12 {f N ■ 75) + V LS (r){L ■ S+) + V LST (r){L ■ S+)(f N • 75) 

+V ALS (r)(L ■§_) + V ALSr (r)(L ■ S„)(t n • r H ) + 0(V 2 ). (3) 

Here S12 = 3(c?i ■ n){(j2 ■ ft) — <7i ■ d% is the tensor operator with ft = f/\f\, S± = (c?i ± c?2)/2 are symmetric 
(+) and antisymmetric (— ) spin operators, L = —irx V is the orbital angular momentum operator, and r/y 
(t~) is isospin operator for N = (p, n) T (S = (S°,5 _ ) T ). We note that the antisymmetric spin-orbit forces 
(Vals and ValSt) do not arise in the NN case because of the identical nature of the nucleon within the 
isospin symmetry. 

According to the above expansion, the wave function should be classified by the total isospin /, the total 
angular momentum J = L + S+ and the parity. A particular spin (isospin) projection is made in terms of 
<7zv-c?h (tjv • th), e.g., for the isospin projection we have p( /=0 ) = (1 — tn -ts)/4 and p( /=1 ) = (3 + rjv7s)/4. 

The equal-time BS wave function for (I,I Z ) = (1, 1) and L = (5- wave) on the lattice is obtained by 



W = 24 S 1? E P ^ (° \p°( r W + k) , (4) 

Pa{x) =S a bc {u a {x)Cj 5 d b {x)) U ca (x), (5) 

^(y) = £abc {u a (y)Cj 5 s b (y)) s cP (y), (6) 



where a and (3 denote the Dirac indices, a, b and c the color indices, and C = 7472 the charge conjugation 
matrix. The summation over R g O is taken for cubic transformation group to project out the 5*-wave statcFH 
The summation over x is to select the state with zero total momentum. Here we take local field operator 
p a (x) and Sjg(y) for the proton and 5°. The wave function and the potential (or equivalently the off-shell 
behavior of the scattering amplitude) depend on the choice of interpolating operators. This is the situation 
common to any field theories. In this paper, we focus exclusively on the local operators as introduced above 
and leave further discussions on the operator dependence to [2S]. We take the upper components of the 

Dirac indices to construct the spin singlet (triplet) channel by P^~ ^ — (02)0/3 (-f > l/3 _1 ' > = ( <j i)q/3)- The BS 
wave function <f>(f) is understood as a probability amplitude to find "nucleon-like" three-quarks located at 
point x + r and "H-like" three-quarks located at point x. Our BS wave function has information not only 
of the elastic amplitude NE — ■> iVS but also of the inelastic amplitudes: The inelastic effects are localized 
in coordinate space at low energies and do not affect the asymptotic behavior of 4>{r) and hence the phase 
shift. 

In the actual simulations, the BS wave function is obtained from the four-point correlator, 

Ppso (x, y, i; to) = (0 \ Pa (x, t)E°p(y, t)J pS o (t ) \ 0) (7) 
= J2 An (° \P^H(V)\ E n) c- £ " (t - to) . (8) 

n 

Here J p s°(to) is a wall source located at t = to, which is defined by J pH o(£o) = Pa/3Pa{to)^{to) with p a (ta) = 
Y,x u x 2 ,s 3 £abc(ua(xi,t )Cj 5 d b (x 2 ,to))u ca (x 3 ,t ) andS^(t ) = T,ff u p 2 ,ff 3 £abc{Ua{yi,to)Cj 5 Sb(y2,to))s c p(y 3 ,to) 
The cigen-encrgy and the eigen-state of the six quark system are denoted by E n and \E n ), respectively, with 
the matrix element A n (to) = (E n \ J p s° (to)\0). For t — to 1, the F pS o and hence the wave function <f> are 
dominated by the lowest energy state. 

The lowest energy state \Eq) created by the wall source J p ~,o(to) contains not only the S'-wave A r S 
component but also the components which can mix by QCD dynamics, such as the D-wave due to the tensor 

3 Due to the periodic boundary condition, this projection cannot remove the higher orbital components with L > 4, which 
however are expected to be small in the ground state. 
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force and AS state due to the rearrangement of the quarks inside baryons. In principle, these components 
can be disentangled by preparing appropriate operator sets for the sink. Study along this line to extract the 
mixing between the S- wave and the D-wave in low energy NN interaction was put forward recently in |28| . 
In the present paper, instead of making such decomposition, we define an effective central potential Vc{r) 
for the (J?- wave component according to our previous works [16|17|18] : 

1 V 2 d>(r) 

Such a potential in the 1 = 1 sector has only the spin dependence as V^ /_1 '(r) = Vo(r) + V a (r) (<tjv • <?h) 
where effects of the tensor force and coupled-channel effects are implicitly taken into account in Vo i(T . 
It is in order here to make some remarks on the potential we have defined. 

(i) The description by the energy-independent non-local potential U (r, f ) is equivalent to that by the 
momentum-dependent local potential, V(r, V). The potential Vc{r) in Eq.© is the leading order 
term of the derivative expansion of V(f, V) and can be determined by the BS wave function at E ~ 0. 
Higher order terms can be extracted successively by the BS wave functions with different E. It was 
recently reported that momentum-dependence of V(f, V) is small at least up to the center of mass 
momentum p = 250 MeV for the NN system J^Hj . 

(ii) The non-local potential U(r,r t ) has one-to-one correspondence to the baryon interpolating operator 
adopted in defining the Bethe-Salpeter amplitude. Different choices of the interpolating operator would 
give different BS wave function and the baryon-baryon potential, although the phase shifts and binding 
energies are unchanged. An advantage of working directly in lattice QCD is that we can unambiguously 
trace the relation between the baryon interpolating operator and the BS wave function or equivalently 
the associated potential. 



3. Lattice setup and hadron masses 



We use the plaquette gauge action and the Wilson fermion action with the gauge coupling )3 = 5.7 on the 
32 3 x 32 lattice. The heatbath algorithm combined with the overrelaxation is used to generate quenched gauge 
configurations. After skipping 3000 sweeps for thermalization, measurement is made for every 200 sweep. The 
Dirichlet (periodic) boundary condition is imposed for quarks in the temporal (spatial) direction. The wall 
source is placed at to = 5 with the Coulomb gauge fixing. These setup are same as our previous calculation 
for the NN potential [16] . 

Masses of pseudo-scalar and vector mesons obtained at hopping parameter (m, k 2 ) are fitted as 

(rop .^ = |(±-l) + |(i_i), v = c+ |(±_±) + |(±_±). (10) 

2 \Kl K C J 2 \K2 K C J 2 \Kl K C J 2 \K 2 K c J 

for Ki = 0.1678,0.1665,0.1640 (i = 1,2). The fit gives k c = 0.16930(1) for the critical hopping parameter. 
From = 135 MeV and m p = 770 MeV, we determine the hopping parameter for physical ud quarks and 
the lattice spacing as K p h ys = 0.16910(1) and a = 0.1416(9) fm (1/a = 1.393(9) GcV), respectively, while 
the hopping parameter for the strange quark mass is given by k s = 0.16432(6) from ijik = 494 MeV. A 
corresponding lattice volume is (4.5fm) 4 , which is large enough to accommodate two baryons. 

In order to check the thresholds of two baryon systems with strangeness S = —2 (AA,A^S,AE and ££), 
we calculate the two-point correlator, C(t;t ) = J2z(®\Ba{x,t)JB a (to)\Q), for the octet baryons (B = 
N, 5, A, E), and Js Q (to) is the wall-source for B. The interpolating fields for A and S + employed in this 
work are 

A Q (x) = e abc {(d a (x)C'j 5 s b (x)) (x)Cj 5 u b (x)) d ca {x) - 2 (u a {x)C~f 5 d b (x)) s ca (x)} , (11) 

S /3 (y) = ~ £ «bc (u a (y)Cj 5 s b (y)) u c p(y). (12) 

Table [1] lists the meson and baryon masses measured at two values of k uc i, 0.1665 and 0.1678, with fixed 
K s = 0.1643. At K u d = 0.1678 corresponding to ~ 368 MeV, 17 exceptional configurations are excluded 
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Fig. 1. (Left) The radial wave function of pS°, in 1 5o (circle) and 3 Si (triangle) channels, obtained at i — to = 6. The Green's 
functions G(r,k 2 ) with r = (r, 0, 0) in 1 So (solid line) and 3 Si (dotted line) are also shown. The inset shows its enlargement. 
(Right) The effective central potential for pS°, in the 1 So (circle) and 3 Si (triangle), obtained from the wave function at time 
slice t — to = 6. The inset shows its enlargement. 

from totally 1300 gauge configurations for the average, while no such exceptional configuration appears 
at K u d = 0.1665 (to t ~ 511 MeV). As seen in Table [TJ the present results for the baryon masses are 
consistent with the experimentally observed ordering of the two-baryon thresholds in the strangeness —2 
sector: i? th (AA) < E th (NE) < E t h(AS) < E t b(ET,). This warrants that AS in the 1=1 channel treated in 
this paper is indeed the lowest energy scattering state. 



4. Results of TVS interaction 



4.1. BS wave function 

The left panel of Fig. [T] shows the wave functions obtained at time slice t — to = 6 with to = 5. The open 
circle (triangle) corresponds to the case of the 1 S'o ( 3 Si) channel. They are normalized to be unity at the 
spatial boundary f = (32/2, 0, 0). For r < 0.7 fm, the wave function is calculated for all possible values of r, 
while, for the outer region, it is calculated only on the x, y, z axis and their nearest neighbors to reduce the 
numerical cost. 

As seen in the Figure, the wave functions are suppressed at short distance and arc slightly enhanced at 
medium distance in both 1 Sq and 3 Si channels. There is a sizable difference of the suppression pattern at 
short distance between and Si; the repulsion is stronger in the channel. 

4.2. central potential 

To derive the AS potential through Eq. ([9]), we need to know the non-relativistic energy E = k 2 / (2/j). It 
has been shown in Ref. [24] that k 2 and thus E can be accurately determined by fitting the wave function 
in the asymptotic region in terms of the lattice Green's function 



K ud k s N con{ m p m K m K * rn^ m p m s o ttia m E+ 

0.1665 0.1643 1000 511.2(6) 861(2) 605.3(5) 904(2) 946(1) 1300(4) 1419(4) 1354(4) 1375(4) 
0.1678 0.1643 1283 368(1) 813(4) 554.0(5) 884(2) 946(1) 1167(7) 1383(6) 1266(6) 1315(6) 
Exp. 135 770 494 892 1019 940 1320 1116 1190 

Table 1 

Hadron masses in the unit of MeV calculated for the hopping parameters K ud and k s , with the number of gauge configurations 
N i 
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Fig. 2. The t dependence of the potential in the 1 So channel (left panel) and the 3 Si channel (right panel) for several values of 
r. The horizontal lines denote the central values at t — to = 6 which arc adopted for the potential in Fig.l (right). The results 
indicated by [6 — 10] are obtained by a constant fit to the data for t — to = 6 — 10 with a single-elimination of the jackknifc 
method. 

G ^ fc2 ) = h E ^p e "' r = {p ] p = z3 j ' (13) 

per y y ' 

which is the solution of (A + k 2 )G(f, k 2 ) = — Jt(r) with 5l{t) being the periodic delta function [12124] . 
Results of the fit in the range (12 < x < 16, < y < 1, z = 0) are shown in the left panel of Fig. [TJ The 
fitting range is determined so that it is outside the range of the interaction (see the right panel of Fig.Q]). The 
non-relativistic energies thus obtained in the ( 3 S'i) channel becomes E = -0.4(2) McV (E = -0.8(2) 
MeV). We have checked that the results of the fit in different ranges 11 < x < 16 and 13 < x < 16 introduce 
systematic errors only less than half of the statistical errors for E. Note that E can be negative for the 
scattering state in a finite box if the interaction is attractive. 

The effective central potentials for pE° system obtained from Eq. © with the wave function <f> and the 
energy E at t — t — 6 are shown in the right panel of Fig. [T]for the 1 S'o and 3 Si channels at m w ~ 368 
MeV. In order to check whether the ground state saturation of the pEP system is achieved in the present 
results, we plot the i-dependence of the potential at distances r = 0.14, 0.20, 0.71, 1.42 and 2.27 fm in Fig.[2j 
The fact that the potential is stable against t for t — to > 6 within errors indicates that the ground state 
saturation is indeed achieved at t — to = 6. This is the reason why we adopted the values at t — to = 6 in 
the right panel of FigQ] Note that a constant fit to data for t — to = 6 — 10 with a single elimination of the 
jackknifc method docs not introduce appreciable change of the final values and the errors of the potential 
as shown in Fig[2j 

In the right panel of Fig.[T] the potential in the iVS system shows a repulsive core at r < 0.5 fm surrounded 
by an attractive well, similar to the NN system |16ll7j . In contrast to the NN case, however, one finds 
that the repulsive core of the pS° potential in the 1 5 f channel is substantially stronger than that in the 
3 S'i channel. Such a large spin dependence is also suggested by the quark cluster model [B]. The relatively 
weak attraction in the medium to long distance region that (0.6 fm < r < 1.2 fm) is similar in both 
and 3 S'i channels. As the energies determined above and scattering lengths given later indicate, the present 
potentials arc weakly attractive on the whole in both spin channels in spite of the repulsive core at short 
distance. Also, the attraction 3 Si channel seems a little stronger than that in the 1 S'o channel. 



4.3. quark mass dependence 

Figure [3] compares the pS° potential at ~ 368 MeV with that at ~ 511 MeV in the 1 S'o channel 
(left) and in the 3 S*i channel (right). At ~ 511 MeV, the potentials are evaluated at t — to = 7. The 
non-relativistic energy in this case is E = —0.19(4) MeV (—0.34(4) MeV) in the 1 S'o ( 3 <Si) channel, where 



G 



> 

X 



m„=368MeV 
05 m„=511MeV 
ffl OPEP(a=0.36) 

% 
P 

'So 




m^368MeV i-^- 
mj=511MeV >— e- 
OPEP(ct=0.36) 




Fig. 3. (Left) The effective central potential for pH° in the x So channel at m-n- ~ 368 MeV (circle) and m n ~ 511 MeV (box). 
The central part of the OPEP (F/{F + D) = 0.36) in Eq. dl4l is also given by solid line. (Right) Same as the left figure, but 
in the 3 Si channel at ~ 368 MeV (triangle) and ~ 511 MeV (diamond). 

the region that (10 < a; < 16, < y < 1, z = 0) is adopted to fit the wave function by the Green's function. 
We have checked that the results of the fit in different ranges 9 < x < 16 and 11 < x < 16 introduce 
systematic errors only less than half of the statistical errors for E. 

The height of the repulsive core increases as the ud quark mass decreases, while the significant difference 
is not seen in the medium to long distances within the error bars. 



4.4. one-pion exchange 

The exchange of a single ir° would conduct the interaction at long distance in the pEP system: 



V* - (1 2a) 9 * NN ^ N ' ^ ( m ^ \ — 

\2m N J 



The pseudo- vector ttNN coupling f^NN and the ttES coupling f^ss are related as / t hh = —f-irNN^ — 2a) 
with the parameter a = F/(F + D) [5]. Also we define j^jv = f-nNN^^- 

The solid lines in Figj3]is the one pion exchange potential (OPEP) obtained from Eq. (fT4]) with m w ~ 368 
MeV, toat ~ 1167 MeV (corresponding to K u d = 0.1678) and the empirical values, a ~ 0.36 [30J and 
g^ NN /(An) ~ 14.0 [I]. Unlike the A^A^ potential in the S- wave, the OPEP in the present case has opposite 
sign between the spin-singlet channel and the spin-triplet channel. Also, the absolute magnitude of OPEP 
is weak due to the factor 1 — 2a. As is seen from FigEl we do not find clear signature of OPEP at long 
distance (r > 1.2 fm) in our potential within statistical errors. On the other hand, there is a clear departure 
from OPEP at medium distance (0.6fm < r < 1.2fm) in both 1 Sq and 3 5i channels. These observations may 
indicate a mechanism of state-independent attraction such as the correlated two pion exchange. 

It should be mentioned here that there is in principle a quenched artifact to the baryon-baryon potentials 
from the flavor singlet hairpin diagram (the ghost exchange) [33j . Its contribution to the central potential has 
a spin-dependent exponential tail Vq (r) oc cr/v ■ <?h exp(— m ff r), which dominates over the Yukawa potential 
at large distances. Its significance can be estimated by comparing the sign and the magnitude of e m ^ r Vc(r) 
in the spin-singlet and spin-triplet channels. Our present data at lightest quark mass = 368 MeV shows 
no evidence of the ghost contribution at large distances within errors. This may indicate the weak coupling 
of the r\ to N and S. 

4.5. scattering length 

Using the energy E = k 2 / (2/i) obtained from the asymptotic behavior of the wave function in Sec l4.21 the 
AS scattering lengths can be deduced from the Liischer's formula |12|24j . 
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Fig. 4. The scattering length for pS° in the 1 Sq (circle) and 3 Si (triangle), obtained by the Liischer's formula from the 
asymptotic behavior of the wave function. The vertical dashed line is the physical point at m n = 0.135 GeV. 



k COt So(k) 



Z 00 (l;<? 2 ) = l/d + O(fc 2 ), 



(15) 



where Zoo(l;g ) with q = is obtained by the analytic continuation of the generalized zeta-function 
■Zoo(s; q 2 ) = -j= Snez 3 ( n2 — <f)~ s defined for Re s > 3/2. The sign of the S*-wave scattering length do is 
defined to be positive for weak attraction. The results arc plotted in Fig0]for two different values of the pion 
masses, ~ 368 MeV and m w ~ 511McV. Both dot (the scattering length in the 3 Si channel) and a us (the 
scattering length in the 1 S'o channel) slightly increase as the quark mass decreases: ao s = 0.05(1) — > 0.10(6) 
fm and dot = 0.09(1) — > 0.21(6) fm. The positive sign indicates that the pE° interaction is attractive on 
the whole in both 1 S'o and 3 Si channels. Due to the possible non- linear dependence of the scattering length 
as a function of the quark mass as is well-known in the case of the NN system [21] (see also [13132) and 
references therein), we are not able to obtain a reliable estimate of the scattering length of the pE° system 
at the physical ud quark mass. In addition, a possible contamination from the ghost exchange may arise at 
small quark masses in the quenched QCD simulations [33] . Therefore, it is necessary in the future to carry 
out the full QCD calculation toward lighter quark masses to predict the precise value of do- 

It would be interesting here to summarize diverse results on the NE scattering length in other approaches. 
The pE° interaction in chiral effective field theory [TIT] predicts weak repulsive scattering lengths of do s ~ — 0.2 
fm and dot ~ —0.02 fm. The phenomenological boson exchange model (e.g., SC97f) [5] gives do s = —0.4 fm 
and dot = 0.030 fm. The quark cluster model (fss2) [7] gives do., = —0.3 fm and dot = 0.2 fm, while QCD 
sum rules gives a$ s = 3.4 ± 1.4 fm and dot = 6.0 ±1.4 fm. 



5. Summary 

We study the pE interaction in the 1 S'o and 3 Si channels through the equal-time Bethe-Salpeter am- 
plitude measured by the quenched lattice QCD simulations on a (4.5 fm) lattice with the quark masses 
corresponding to m x /m p ~ 0.45 and 0.59. We adopt specific choice of the baryon interpolating operators 
as given in Eqs. §5§ and ^ and extract the potential associated with these operators. The effective central 
potential deduced from the Bethe-Salpeter wave function has repulsive core at short distance surrounded 
by attractive well at the medium and long distances. The scattering lengths determined by fitting the 
asymptotic wave function with the solution of the Helmholtz equation indicates that the pE° interaction is 
attractive on the whole in both 1 5'o and 3 Si channels. There is a slight tendency that the 3 Si interaction is 
more attractive than that in the ^o channel. 

To reduce the uncertainties due to the lattice discretization at short distances, we need to carry out 
the simulations with smaller lattice spacing and/or with the improved lattice action. Systematic studies in 
various other channels such as AiV, EiV and AA are not only interesting by themselves but also important for 
studying the structure of hypcrnuclci and the interior of neutron stars. Also, (2+l)-flavor QCD simulations 
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with realistic quark masses will be the ultimate goal to make qualitative predictions of the NN and YN 
interactions: Studies along this line with the PACS-CS gauge configurations [34] are now under way. 
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